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ABSTRACT 



Aims. We study the impact of galaxy-galaxy interactions on the far-infrared properties of galaxies and its evolution at < z < 1.2. 
I ■ Methods. Using the high-z galaxies in the fields of Great Observatories Origins Deep Survey (GOODS) observed by Herschel Space 

P ' Observatory in the framework of the GOODS-Herschel key program and the local IRAS or AKARI-selected galaxies in the field of 

Sloan Digital Sky Survey Data Release 7, we investigate the dependence of galaxy properties on the morphology of and the distance 
C/2 , to the nearest neighbor galaxy. 

Results. We find that the star formation rates (SFRs) and the specific SFRs (SSFRs) of galaxies, on average, depend on the morphology 
of and the distance to the nearest neighbor galaxy in this redshift range. When a late-type galaxy has a close neighbor galaxy, the SFR 
and the SSFR increase as it approaches a late-type neighbor, which is supported by Kolmogorov-Smimov (K-S) and Monte Carlo 
(MC) tests with a significance level of >99%. However, the SFR and the SSFR decrease or do not change much as it approaches an 
early-type neighbor The bifurcations of SFRs and SSFRs depending on the neighbor's morphology seem to occur at R„ a 0.5rvi,.nei 
(virial radius of the neighbor), which is supported by K-S and MC tests with a significance level of >98%. For all redshift bins, the 
SSFRs of late-type galaxies interacting with late-type neighbors are increased by factors of about 1.8±0.7 and 4.0±1.2 compared to 



both local and high-z late-type galaxies strongly interacting with late-type neighbors (i.e. R„ < 0. lrvi|.„ei) appears to be higher than that 
0^ ' of non-interacting galaxies with a significance level of 96 - 99%. However, the dust temperature of local late-type galaxies strongly 

' interacting with early-type neighbors seems to be lower than or similar to that of non-interacting galaxies. 

Conclusions. Our results suggest that galaxy-galaxy interactions and mergers have been strongly affecting the SFR and the dust 
properties of star-forming galaxies over at least 8 billion years. 

Key words, galaxies: active - galaxies: evolution - galaxies: formation - galaxies: interactions - galaxies: starburst - infrared:galaxies 
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1 . Introduction on the distance to and the morphology of the nearest neighbor 

galaxy. When a galaxy is located within the virial radius of its 
In the hierarchical picture of galaxy formation, galaxies are nearest neighbor, its morphology tends to be the same as that of 
formed and grow through continuous interactions and mergers j^e neighbor This indicates an important role of hydrodynamical 
with other galaxies. These galaxy-galaxy interactions and merg- interactions with neighbors within the virial radius. This mor- 
ers are expected to strongly affect galaxy properties such as mor- phological confor mity was also found between host and their 
phology, luminosity, structure parameters, star formation rate satellites galaxies (lAnn et al .''2008: Wang et al. 2010b). and be- 
(SFR), or dust properties (see Struck 2006 for a review). j^een galaxies in galaxy clusters (Park & Hwang 2009 ). 

Since the merger between spiral galaxi es is known t o 
form elliptical galaxies as first suggested by iToomrd (Il977h . Fo r the star formation activity (SFA), iLars on & TinslevI 
there is growing evid ence for a change in galaxy morphology (Il978h first noted that morphologically normal and peculiar 
(e.g. Park et al. 2008) and galaxy structure with merg er (e.g., galaxies show very different optical color distributions, which 
Nikolic et al. 2004; Patton et al. 2005; Hernandez-Tol edo et al.l suggests that the SFA of peculiar galaxies is enhanced com- 
2005[ iPark&Chol [20091) . For example, iPark et air(l2008l) pared to those of normal galaxies. Many studies have exten- 



showed that galaxy morphology and luminosity strongly depend sively investigated this enhancement of SFA in paired galax- 
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ies (e . g.. ICondon et al.lll982l: iKeel et alJll985t iKennicutt et al] 
1987"; 'Barton et al.' '200d^, 'Lambas et al.l_|2003|; iNikolic et al.' 
2004 ; Alonsoetal. 2004; De Propris et a! ] 
20061: IWoods & Geller 2007; Li et al. 2008; 




Knapen & James 20 09; ,Perez et al., ,ZUU^; ,FarK & Choii 
Hwang etal. 2010a; " iXu et al.1 l2010riDarg et alj 120101 



m 

200l 'Geller et al. 

Ellison et al. 
2009: .Park & Choi. 



Patton et al.ll201 li) . In summary, the SFA of galaxies seems to 



strongly depend on the morphology and the mutual separation of 
galaxies in pairs. For example, the SFR of a galaxy is increased 
when the nearest neighbor has a late morphological type, but is 



decreased or remains the same when the neighbor is an early- 
type galaxy (jPark & Choil 120091; iHwang et al]|2010at IXu et alj 
12010 ). There are two characteristic pair-separation scales where 
the SFA abruptly changes: the virial radius (rvir,nei, to be de- 
fined in g2.4b of the nearest neighbor galaxy where the effects of 
galaxy interaction emerge (roughly a few hundred kpc for bright 
galaxies ) and ~0.05 ryime i where the galaxies in pairs start to 
merge ( Park&Choill2009l) . However, some other studies found 
no significant difference i n galaxy properties between isolated 
and paired galaxies (e.g., ^ Yee & EUingsonl [19951: iPatton et"al] 
[l997; Bergvall et al. 2003). For high-z galaxies, the evolution 
of the fraction of interacting galaxies has been st udied exten- 
sively and is found to increase with redshift (e.g ., IZepf & Kool 
Il989t iKartaltepe et all 120071; [Bridge et alj|201o[) . but there are 
few studies focusing on the enhancement of SFA in paired galax- 
ies (e.g., Lin et al. 2007; de Ravel et al. 2009; Woods et al. 2010; 
IWong et al.ll20ni) . 

For dust-enshrouded systems such as luminous infrared 
galaxies (LIRGs), one might have expected that not only the 
SFR but also the dust properties such as dust temperature or 
dust mass change as a function of pair separation. Interestingly, 
iTelesco et al.l (Il988h found the highest dust color temperature 
(i.e., the flux density ratio between IRAS 60 and 100 fim) for 
galaxies in pairs w ith the smallest separation. O n the other hand, 
iLutz et all (Il998h andlRigopoul ouet alJ 11999) did not find any 
correlation between the line-to-continuum ratio for the 7.7 fj.m 
PAH emission feature (i.e. SF vs. AGN activity) and the nu- 
clear separation of the interacting components of ultraluminous 
infrared galaxies (ULIRGs). Klaas et al. (2001) also searched 
for the differences in the mid-infrared (MIR), far-infrared (FIR) 
and submm luminosity ratios of interacting and non-interacting 
ULIRGs, but did not find any noticeable difference depending 
on the merging stage. HoweverJ Xilouris et al.l (12004') found evi- 
dence for an increase in the 100- to 850-/zm flux density ratio (a 
proxy for the mass fraction of the warm and cold dust) with the 
merging sequence. They argue that the contr adictory resii l ts may 
be caused by the different sample selection. iKlaas et al] (1200 lb 
did not order their interacting ULIRG systems on the basis of the 
pair separation. Their samples span a small range in bolometric 
luminosity, therefore, they might have missed some interacting 
systems in the ver y early stage of merging. It is also noted that 
ISmith et al.1 (l2007h found no significant change in Spitzer MIR 
colors with pair separation in their sample of tidally distorted 
premerger interacting galaxy pairs. This may be because their 
sample was selected to be in the early stage of interactions. 

Thanks to the advent of the Herschel Space Observatory 

dPilbr att et al. 2010) with its very wide wavelength coverage 
(70 - 500/zm), we are now able to have both the 'Wien' 
and 'Rayleigh-Jeans' sides of the FIR spectral energy distri- 
butions (SEDs) of high-z galaxies, and to have robu st SFR 
and dust temperature measurements for them (e.g., Elbaz et al.' 
[2010, 2011; Hwang etal. 2010b; Magdis et al. 2010ab). In 
this paper, we investigate the impact of galaxy-galaxy interac- 
tions on galaxy properties focusing on FIR properties of high- 



z galaxies at 0.4 ^ z ^ 1.2 by taking advantage of the 
wide wavelength coverage of the Photodetector Array Camera 
(PACS; Poglitsch et al. 2010) and Spe ctral and Photometric 
Imaging Receiver (SPIRE; Grif fin et aiDl2010t) instruments on- 
board Herschel in the fields of Great Observatories Origins Deep 
Survey (GOODS; [Dickinson et alj2003l:lGiavaUsco et alj20"04l) . 
To compare these results with those of local galaxies, we also 
use a sample of galaxies at z < 0.1 that were covered by IRAS 
anAAKARI satellites (Murakami et al. 2007) in the field of Sloan 
Digital Sky Survey (SDSS, ,York et al.. ,2000ft . The IRAS and 
AKARI all-sky survey data contain flux density measurements 
at 12-160 jum, which can probe both the 'Wien' and 'Rayleigh- 
Jeans' sides of the peak of IR SEDs of local galaxies in a way 
similar to that of Herschel for high-z galaxies. Section [2| de- 
scribes the data used in this study, and the dependence of galaxy 
properties on the distance to and the morphology of the near- 
est neighbor galaxies is given in |3] Discussion and conclu- 
sions are given in ^and ^ respectively. Throughout, we adopt 
h - 0.7 and a flat ACDM cosmology with density parameters 
Qa,o = 0.73 and Q„, o = 0.27. 



2. Data 

2.1. GOODS Sample 

We used a spectroscopic sample of galaxies in GOODS, which is 
a deep multiwavelength survey covering two carefully selected 
regions including the Hubble Deep Field North (HDF-N) and the 
Chandra Deep Field South (CDF-S). Hereafter, the two GOODS 
fields centered on HDF-N and CDF-S are called GOODS-N and 
GOODS-S, respectively. The combined area of the two fields is 
approximately 320 arcmin^. 

The GOODS fields were observed by Herschel in the 
GOODS-He rschel key program (Elbaz etal. 2011). The fufl 
10' X 16' GOODS-N field was observed with both PACS (100 
and 160 pm) and SPIRE (250, 350, and 500 pm). A smaller re- 
gion within the GOODS-S field (x 10' x 10') was observed with 
PACS only. We supplemented these data with public SPIRE im- 
ages covering the full 10' x 16' GOODS-S field, which was orig- 
inally taken in the HerMES key program (Oliver et al. in prep.). 
Source extraction on these PACS and SPIRE images was per- 
formed at the prior positions of Spitzer 2 4 /im-s elected sources, 
and details are described in lElbaz et aP (1201 ll) . Note that this 
extraction method with 24 ^im-s elected sources reco vers more 
than 99% of Herschel sources (iMagdis et al.l l20lH) . We used 
flux densities in PACS bands down to 3cr limits of 1.1 and 
2.7 mJy (0.8 and 2.4 mJy) at 100 and 160 pm in GOODS-N 
(GOODS-S), respectively. SPIRE measurements are used down 
to So- limits of 5.7, 7.2, and 9.0 mJy at 250, 3 50 and 500 pm, 
respectively (see Table 1 in lElbaz et al.ll20ll1 for more details 
about the noise properties). By combining Herschel data with 
the existing multi-wavelength data, we made a band-merged cat- 
alog of GOODS galaxies using the p hotometric d ata at HST 
ACS BViz, CFHT/WIRCam JK (Norfli;IWang et al.l[2010 a) and 



201^, Spitzer IRAC 



VLT/ISAAC JHK (South; Ret zlaff et al.l 
3.6, 4.5, 5.8, 8 pm and MIPS 24 and 70 pm, Herschel PACS 
100 and 160 pm, and Herschel SPIRE 250, 350 and 500 pm. 

Among the sources in the band-merged catalog, we 
used 3630 and 3181 galaxies whose spectroscopic red- 
shi fts are reliable over the whole fi e lds of GOODS- 
N (Cohen et al.' '2000; 'Cowie et al.' '2004'; 'Wirth et al.' '2004 
Reddyetal. 2006; BargeretjiL .2008; Cooper etal. 20ll 
Stern et al. in prep) and GOODS-S iSzokoly et al. 2004 
—i\ 



iLe Fevre et al.ll2004 [MignoU et alJl2005UVanzeUa et aL,2005 
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Fig. 1. (a-b) Dust temperature (Tdust), (c-d) IR luminosity, (e-f) stellar mass, and (g-h) evolution-corrected, rest frame r-band absolute 
magnitude vs. spectroscopic redshift for galaxies in (left) SDSS and (right) GOODS. Red dots indicate FIR-selected galaxies. Black 
dots denote galaxies without FIR detection in the spectroscopic sample of galaxies (only 3 % of SDSS galaxies in the total sample 
are shown). Open blue circles denote galaxies with Tdust measurements. Blue solid lines in (g-h) define the volume limited samples 
of SDSS galaxies at 0.005 < z < 0.0656 and of GOODS galaxies at 0.4 < z < 1.2, respectively. The bottom curve in (g) indicates 
the a pparent magnitude limit (nir - 17.77) for the main galaxy sample in SDSS using the mean iiT-correction relation given by Eq. 
(2) of lChoi etal.l (l2007h . 



20061 I2008t iRayikumaretalJ 1200 

Kurk et alj 120091: iBalestra et alj 1201 Ot [Silverman etall 



Popesso et al. 



2009; 



201C; 



Xia et alJ l201 ih . respectively. The rest frame r-band abso 



lute magnitude Mr of galaxies was computed based on the 
ACS plus near-inf rared (NIR) photometry with TiT-corrections 
dBlanton & RoweiSi,2007,) . It was computed in fixed bandpasses, 
shifted to z - 0.1 to be compared with local SDSS galax- 
ies. The l. l(z - .1) term was added to M,- for the evolution 
correction (IWolfet al. 2003). We calculated stellar masses us- 
ing t he photometric data upto IRAC 4.5 yum with a Salpeter 
IMF (ISalpeteiilI955l) . We use d the code Z -PEG with a galaxy 
evolution model PEGASE.2 (iFioc & Rocc a-Volmerange 1999; 
iLe Borgne & Rocca-Volmerangell2002l : see lElbaz etalj|2011l for 
details). 

We adopted the galaxy morph ology information [earl y types 
(E/SO) and late types (S/Iit)] from lHwang & ParkI (|2009|) that is 
based on the visual inspection of ACS BViz images. We per- 
formed additional visual classification for the galaxies in the 
volume-limited sample of gal axies shown in the right panels of 
Fig.[T]that are not included in lHwang & Pa rk (2009). 

We computed the IR luminosity (Lir) for 903 and 828 galax- 
ies with spectroscopic redshifts in GOODS-N and -S, respec- 
tively, detected in at least one out of the five PACS/SPIRE bands. 
We fit the flux densities at /Irest > 30yum by allowing the normal- 
ization of the SED templates of Chary & Elbaz (2001, CEOl) to 
vary and choosing the one that minimizes the values. When 
there are two or less data points to fit (i.e. < 2), we fit the flux 



densities without allowing the normalization of the templates 
(i.e. standard CEOl technique). There are some cases with only 
one FIR band used for Lir measurement, but the IR luminosities 
extrapolated from a single passband were found to agree very 
well with those measured with all FIR bands with an average un- 
certainty of -30% (see Fig. 23 of lElbaz et aDl201 ih . Therefore, 
this does not introduce any bias in our results. To determine the 
dust temperature, we fit the observational data with a modified 
black body (MBB) model by fixing the emissivity parameter to 
y6 = 1 .5. We require at least one flux measurement at each ' Wien' 
and 'Rayleigh-Jeans' s ide of the FIR pe ak [i.e., at least two mea- 
surements in total; see lHwang et al.ll20 10b for detailed selection 
criteria]. In the result, we have 284 and 104 galaxies with Tdust 
measurements in GOODS-N and -S, respectively. 



2.2. Sloan Digital Sky Survey Sample 

For local galaxies, we used a spectroscopic sample of galax- 
ies in the SDSS Data Release 7 ( Abazajian et al. 2009, SDSS 
DR7) complemented by a photometric sample of SDSS galax- 
ies whose redshift information is not available in the SD SS 
database, but available in the literature dHwang et al.ll2010ah . In 
addition, we used the FIR data for these galaxies compiled in 
Hwang et al. (2010b) by cross-c oiTel a ting I RAS Faint Sources 
Catalog - Version 2 (Moshir etall |1992[) a nd AKARI/Fai- 
Infrared Surveyor (FIS; iKawada et alJ ~ 20071) all-sky survey 
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Bright Source Catalogue (BSCQ) ver 1.0 with the SDSS sam- 
ples. 

The r-band absolute magnitude Mr was also computed 
in fixed bandpasses, shifted to z = .1, using Galactic 
reddening correction ( Schlegel et al.l [T998h and /T-corrections 
(Blariton & RoweiJ 2007 ). The evolution correction given by 
iTegmark et all (12004 . £(7) = 1.6(z - 0.1), is also apphed. Note 
that the amount of evolution coiTection is different between 
GOODS and SDSS galaxies. These values are taken from the 
redshift evolution of the characteristic luminosity in the lumi- 
nosity function of local and high-z galaxies separately, so we 
kept different values. Change of these values does not affect our 
conclusions. 

By adopting a method similar to the one applied to GOODS 
galaxies, we computed the IR luminosity of 14444 galaxies 
(among the total sample of 926,748 SDSS galaxies) whose IRAS 
60 fim or AKARI 90 jum flux densities are reliabl^ using the 
CEOl SED templates by allowing normalization of the tem- 
plates. We fit the flux densities without allowing the normaliza- 
tion of the templates if there are two or less data points to fit. 
To determine the dust temperatures, we again fit the observa- 
tional data with a modified black body (MBB) model by fixing 
the emissivity parameter to /3 - 1.5 only for 238 galaxies de- 
tected at AKARI 140 or 160 /vm so that we can have flux den- 
sity measurements longwards of the FIR peak as well as the one 
shortwards o f the peak in a sim ilar way to GOODS-Herschel 
galaxies (see lHwang et aljr2010bl for detailed selection criteria). 

The stellar mass estimates were obtained from MPA/JHU 
DR7 value-added galaxy catalog (VAGC), which are based 
on the fit of SDSS five -band photometry dKauffmann et alj 
l2003H Gallazzi et al."2005V We convert these estimates based on 
Kroupa IMF (Kroupa 2001 ) to a Sal peter IMF by dividing them 
by a factor of 0.7 ( Elbaz et allllool . 

We adopted the galaxy morphology information from the 
Korea Institiite for Advanced Sfii dy (KIAS) DR7 VAGCQ 
dPark & Choi"2005 UChoi etal] |201Cl). We performed additional 
visual classification for the galaxies in the SDSS database that 
are not included in KIAS DR7 VAGC. 



2.3. Comparison of GOODS and SDSS galaxies 

In the right panels of Fig.[T] we plot several physical parameters 
of the GOODS galaxies as a function of redshift, and define a 
volume limited sample to be analyzed (-19.5 > M,- with 0.4 < 
z < 1.2). Similarly, we plot the SDSS galaxies in the left panels 
of Fig. [T] and define a volume limited sample of galaxies with 
-19.5 > M,. and 0.005 < z < 0.0656. The comoving volume for 
this SDSS sample is ~ 2. Ox 10^ Mpc^, which is much larger than 
that for the GOODS galaxies (-2.4x10^ Mpc^ for GOODS-N at 
0.4 < z < 1.2). Note that the FIR detection Hmits for local SDSS 
and high-z GOODS galaxies are not the same even if we fix the 
mass and luminosity of galaxies in ^ However, since we are 
interested in the difference of SEA depending on the morphology 
of and the distance to the nearest neighbor galaxy in a given 
redshift range, the different FIR detection limits between local 
and high-z galaxies do not affect our conclusions. 




-'SO.IRAS, 



Fig. 2. Correlation between dust temperature (Tdust) and flux 
density ratios (56o,iras/5 ioo.iras) for local SDSS galaxies. The 
dashed line represents the best-fit with an ordinary least-squares 
bisector method. 

In Fig. [2] we plot the dust temperature versus FIR flux 
density ratio (^eo.iRAs/'S ioo.iras) for SDSS galaxies, which 
shows a good corre lation between the two quantities (see also 
IChanial et al]l2007|K From th e fit with an ordinary least-squares 
bisector method dlsobe et al. 1990), we derive a transformation 
relation between the two quantities. 



' http://www.ir.isas.jaxa.jp/AKARI/Observation/PSC/ 
Public/RN/AKARI-FlS JSC.Vl_RN.pdf 

^ Flux quality flags are either 'high' or 'moderate' for IRAS sources 
and 'high' for A A^/?/ sources. 

^ http://www.mpa-garching.mpg.de/SDSS/DR7/ 
http://astro.kias.re.kiyvagc/dr7/ 



Tdust (K) = (43.0 + 0.3) + (37.0 + 1.5)log(5 60,IRAs/'5 100,IRAS)-(1) 

Since the number of SDSS galaxies with S(,qjras/S ioo.iras 
is larger than that with Tdust measurements, we transform the FIR 
flux density ratio into the dust temperature using this equation to 
increase the statistics (see Fig.|6]l. 



2.4. Nearest neigiibor galaxy 

To investigate the effects of interactions with the nearest neigh- 
bor galaxy, we first identified the nearest neighbor of a target 
galaxy that is the closest to the target galaxy on the projected sky 
and that satisfies the conditions of magnitude and relative veloc- 
ity. We searched for the nearest neighbor galaxy among galax- 
ies that have magnitudes brighter than Mr - /Wr.target + 0.5 and 
have relative velocities less than Au — |yneighbors ~ ftargetl = 660 
km s"' for early-type target galaxies and less than Av - 440 
km s"' for late-type target galaxies. Because the redshift uncer- 
tainties are larger for the GOODS galaxies, we use velocity dif- 
ference limits that are 10% larger than those that we have al- 
lowed for SDSS galaxies in previous studies (Av - 600 and 
400 km s ' for early- and late-t ype target SDSS galaxies, re- 
spectively, as seen in Fig. 1 of Park et al.ll2008) . The use of dif- 
ferent values for the relative velocity condition depending on 
galaxy morphology is supported by dififerent velocity distribu- 
tions of neighboring g alaxies dependi ng on galaxy morphology, 
as sho wn in Fig. 1 of lParket al.l(l2008h an d Fig. 2 of|Pa rk & CM 
(l2009h for SDSS galaxies and Fig. 3 of iHwang &Parh (i2009l) 
for GOODS galaxies. Since we used volume-Umited samples 
of galaxies with Mr < -19.5, we restrict our analysis to tar- 
get galaxies brighter than /Wr,target - -20 so that their neighbors 
(Mr^nei ^ A^r.target + 0.5) are Searched within the volume-limited 
samples. 
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The virial radius of a galaxy within which the mean mass 
density is 200 times the critical density of the universe (pc), is 
calculated by 



r™ = (3rL/47r)i/^(200p,) 



-1/3 



(2) 



where L is the galaxy luminosity, and y the mass-to-light ra- 
tio. Here, the mass associated with a galaxy plus dark halo 
system is assumed to be proportional to the r-band luminos- 
ity of the galaxy. We assume that the mass-to-light ratio of 
early-type galaxies is on average twice as large as that of late- 
type galaxies at the same absol ute magnitu de M^, which means 
r(earlv)= 2y(late ) [see §2.5 of iPark & Ch oi (2009 ) and §2.4 of 
[Park et al.' (2008) for SDSS galaxies, and §2.3 of lHwang & ParkI 
(i2009.) for GOODS galaxies]. The critical density of the uni- 
verse Pc is a function of redshift z [i.e. pc - 3//^(z)/(87rG)] and 
^miz) = Pb(z)lpciz) = p(l + zf Ipc(z), where ph and p are the 
mean matter densities in proper and comoving spaces, respec- 
tivel}0. The Hubble parameter at z is //^(z) - //o[Q„,,o(l + zf + 

^^/t,o(l + z)^ + fiA.o], where Q,„_o, ^kf.h and Qa.o are the dimen - 
sionless density parameters at the present epoch (|Peebleslll993h . 
Then, the virial radius of a galaxy at redshift z in proper space 
can be rewritten by 

r™(z) = [3rLO„,,o/(800;rp)/{a„,o(l+z)'+Q^,o(l+z)'+QA,o)]'^'.(3) 

The value of mean density of the universe for local SDSS 
galaxies, p = (0.0223 +0.0005)(rL)-2o(/j"'Mpc)-^ was adopted 
where {yL) 20 is the mass of a late-type galaxy with M,. = -20 
(lParketal.11^08) . 

For high-z GOODS galaxies, we computed the mean mass 
density p using the galaxies at z = 0.4 - 1.2 with various ab- 
solute magnitude limits varying from - -16 to -20. We 
found that the mean mass density appears to converge when the 
magnitude cut is fainter than M,- = -17.5, which means that 
the contribution of faint galaxies is not significant because of 
their small masses. In this calculation, each galaxy is weighted 
by the inverse of completeness according to its apparent magni- 
tude and color (see Fig. 1 of lHwang & Pa rk 2009). We obtain 
p = 0.017 and 0.013 (rL)-20 (Mpc'^) for GOODS-N and -S, re- 
spectively, where (yL)-.20 is the mass of a late-type galaxy with 
Mr = -20. According to our formula the virial radii of galaxies 
with Mr - -20 and -21 are 300 and 400 /i"' kpc for early types, 
and 240 and 320 /i"' kpc for late types, respectively. 

The spectroscopic completeness can affect the identification 
of the genuine nearest neighbor, and then the nearest neighbor 
can be seriously misidentified if the completeness is very low. 
Our previous Monte Carlo experiment shows that the fraction of 
the misidentified nearest neig hbor reaches about 50 % when the 
sample completeness is 50% dHwang & Parkll2b09l) . Therefore, 
it is necessary to have survey data with high completeness in 
order not to miss the genuine nearest neighbor Up to now, 
GOODS has the highest spectroscopic completeness (71 - 86% 
for GOODS galaxies at < 23) among several large, deep- 
field surveys with HST images to our knowledge. Therefore, it 
is the best survey data for our analysis, but it should be noted 
that our results could be weakened by this incompleteness. The 
complete ness depends on the ap parent magnitude and color (see 
Fig. 1 of iHwang & ParkI l2009h and also depends on the dis- 
tance between galaxies due to the difficulty in observing galaxies 
close to each other using multiobject spectrograph (MOS). We 
checked the completeness as a function of the projected distance 



to the target galaxy, and found that it does not change with the 
projected distance. It might be because we combined spectro- 
scopic data from numerous references, therefore, the difficulty 
in observing galaxies with small separation using MOS is sig- 
nificantly reduced. Similarly, the redshift information of some 
SDSS galaxies missed by th e SDSS database wa s complemented 
by the data in the literature (iHwang et al.ll20ldah . so there is also 
no bias for local galaxies. 

Note that the nominal FWHMs of the point spread function 
(PSF) are 37"(65 pm), 39"(90 urn), 58"(14 pm), and 61 "(160 



pm) for AKARI bands (iKawada et al 



2007h. 



and 1 .44'(6 pm) 

and 2.94'(100 pm) for IRAS bands (ISanders etalj 12003). The 
corresponding angular size of one viral radius for typical late- 
and early-type galaxies with M, = -21 and z = 0.035, is 11.3' 
and 14.3', respectively. Therefore, though the galaxies in pairs 
are selected from galaxy catalogs in optical bands with high spa- 
tial resolutions, the measured SFRs of galaxies in close pairs 
cannot be clearly assigned to the individual galaxies due to the 
poor spatial resolution of FIR data. For example, local galaxies 
with Mr - -21 and z - 0.035 having neighbors at < 0.06rvir 
(< 0.12rvi,) are not resolved with AKARI 90 pm (IRAS 60 pm), 
because the pair separation is smaller than the FWHM at each 
band. Thus the measured SFRs for these galaxies can indeed 
indicate SFRs of the whole interacting systems. However, our 
results on the increased SFRs (or SSFRs) due to galaxy-galaxy 
interactions to be seen in ^ are not strongly affected by this 
effect because the increased SFRs (or SSFRs) are found to be 
much larger than a factor of two that could be simply due to this 
blending problem of two galaxies in one FIR beam. 

For high-z galaxies, the PSF FWHMs are 6.0" (Spitzer 24 
pm), 6.7" {Herschel 100 yum), 11.0" (160 yum), 18.1" (250 jum), 
24.9" (350 pm), and 36.6" (500 pm). The angular size of one 
virial radius for typical late- and early-type galaxies with Mr - 
-21 would be, respectively, 70.8" and 88.6" at z = 0.6, and 
59.1" and 73.5" at z = 1.0. Similar to the case of local galaxies, 
high-z galaxies having neighbors at < 0.07rvir (< 0.09rvii) are 
not resolved with Spitzer 24 pm (source extraction on Herschel 
images was performed at the prior positions of Spitzer 24 yum- 
selected sources: Elbaz et aL.201 1.) . However, this has no effect 
on our conclusions. 



^ We correct a typo for Q.^{z) in IHwang & ParU ( I2009I) . but all the 
related values in that paper are correct. 



3. Results 

3.1. Change in SSFRs as a function of pair separation 

In Fig. [3] we plot several FIR properties of local (left) and high-z 
(right) late-type galaxies as a function of the projected distance 
to the nearest neighbor galaxy normalized by the virial radius of 
the neighbor (/?n/'"vir.nei)- Note that we only plot IRAS orAKARI- 
selected (left) and Herschel-selected (right) late-type galaxies, 
but their nearest neighbors are selected among the spectroscopic 
samples of SDSS and GOODS galaxies regardless of FIR de- 
tections. To remain complete in terms of neighbor galaxies and 
remove the effect of mass on the SFA, we restrict our analysis 
to the late-type galaxies with -20 > Mr and 10^° < Mjtar (Mq) 
< 10". Therefore, even if we use the physical distance for the 
pair separation without using a normalization by the virial ra- 
dius, the observed trends in this figure do not change. 

For high-z samples, we only plot GOODS galaxies that have 
at most one neighbor within 6" (the full width half maximum, 
FWHM, of the Spitzer beam at 24 pm) with ^24 > 50% of the 
central 24 pm source (i.e. A^nei,24/jm ^ !)■ This criterion is intro- 
duced to reduce the contamination of neighboring sources, but to 
allow for the possibility of having one close neighbor in order to 
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Fig. 3. (a-b) SFR, (c-d) M^tar, (e-f) specific SFR for late-type galaxies with -20 > M,- and 10'° < M,s,a,- (Mq) < 10" as a function 
of the projected distance to the nearest neighbor galaxy (/?n/'"vir,nei) in (left) SDSS (0.005 < z < 0.0656) and (right) GOODS 
(0.4 < Zspec < 1.2). Red and blue dots are late-type galaxies having early- and late-type nearest neighbor galaxies, respectively. 
Large circle (early-type neighbor case) and square (late-type neighbor case) are median values of each physical parameter and 
of /?n/rvir,nei at each distance bin. The errorbars represent 68% (Icr) confidence intervals that are determined by the bootstrap 
resampling method. 



study the effects of galaxy proximity. In practice, this criterion 
results in the removal of only 4 galaxies among 330 GOODS 
galaxies in Fig. [3] which has no effect on our conclusions. 

Assuming a Salpeter IMF, we converted the I R luminosity 
into SFRiR using the relation in iKennicut 3 (Il998h : SFRiR (M^ 
yr-') = 1.72 x 10"'°Lir(Lo). Panels (a-b) show that SFRir de- 
pends on the projected distance to the nearest neighbor as well 
as on the neighbor's morphology. Large circle (early-type neigh- 
bor case) and square (late-type neighbor case) are median values 
of each physical parameter and of rvir,nei at each distance bin. 
It is seen that SFRir changes with neither the distance to the 
neighbor nor the neighbor's morphology when a galaxy is lo- 
cated farther than 0.5 rvir,nei (virial radius of neighbor galaxy). 
On the other hand, when a galaxy is located at /?„ ^ 0.5rvii,nei, 
the SFRir increases as the target, late-type galaxy approaches a 
late-type neighbor, but decreases or remains constant as it ap- 
proaches an early-type neighbor. It is important to note that the 
bifurcation of SFRir depending on the neighbor's morphology 

occurs ati^n ~ 0.5rvir,nei- 

Since the SF R is strongly correlated with the mass 
(Brinchmann et al 2j2004t lElbaz et al.1 l2007t iDaddi et al.1 l2007t 
[Magdis et al.l l2010ch . it is important to check the effects of 
galaxy mass on the change of SFR. Thus we plot the distribu- 
tions of stellar masses and specific SFRs (SFRiR/Mstai, SSFRs) 
in panels (c-f). Panels (c-d) show no significant change of masses 
as a function of R^, but panels (e-f) suggests that a differ- 
ence in SSFRs between early- and late-type neighbor cases is 
still prominent at < 0.5rvii,„ei. We also checked the red- 
shift distributions as a function of (not shown), which shows 
again no significant dependence on /?„ and the morphology of 
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Fig. 4. Distribution of SSFRs (top) and Tdust (bottom) for SDSS 
(left) and GOODS (right) galaxies. Non-interacting galaxies (i.e. 

> '"vir.nei in Fig. O with early- and late-type neighbors are 
shown by solid and dotted histograms, respectively. Interacting 
galaxies (i.e. < 0.5rvii,nei) with early- and late-type neigh- 
bors are denoted by hatched histograms with orientation of 45° 
(// with red color) and of 315° (\\ with blue color) relative to 
horizontal, respectively. 
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Fig. 5. SSFRs of interacting galaxies normalized by SSFRs of non-interacting galaxies (i.e. R„ > ryir.nei in Fig. [3]) as a function of 
redshift for a) early-type and b) late-type neighbor case. Filled and open circles are median values for strongly- (i.e. ^ 0.1rvii,nei), 
and weakly-interacting (i.e. /?„ < O.Srvimei) galaxies, respectively, at each redshift bin (0.005 < z < 0.0656, 0.4 < z < 0.8, and 
0.8<z<1.2). 



neighbor. We used a Kolmogorov-Smirnov (K-S) test to deter- 
mine whether the SSFR distributions of interacting and non- 
interacting galaxies with late-type neighbors are drawn from 
the same distribution. We tested two cases of interacting galax- 
ies with late-type neighbors : 1) strongly-interacting (Rn < 
0.1?"vir,nei) galaxics vcrsus non-interacting galaxies (R^ > ryir.nei), 
and 2) relatively weakly-interacting (/?„ ^ 0.5rvii,nei) galaxies 
versus non-interacting galaxies. The hypothesis that the two dis- 
tributions are extracted from the same parent population can 
be rejected at the confidence level of >99% for both cases in 
the local universe. In the high-z universe, the hypothesis can be 
rejected at the confidence level of >99 and 90% for strongly- 
and weakly-interacting galaxies, respectively. To help the under- 
standing of the difference in SSFRs for several subsamples, we 
plot, in Fig.m the distribution of SSFRs for interacting and non- 
interacting galaxies in SDSS and GOODS. 

The statistical significance of the increased SSFRs of in- 
teracting galaxies with late-type neighbors compared to non- 
interacting galaxies, is also tested by Monte Carlo (MC) test. We 
constructed two subsamples (with the same number of galax- 
ies as in the actual samples) by randomly drawing the SSFRs 
from the whole galaxy sample and computed the median of 
each subsample. The resulting two subsamples will have the 
same medians on average. These random subsamples will tell 
us about whether or not the difference in SSFRs between inter- 
acting and non-interacting galaxies with late-type neighbors is 
statistically significant. We constructed 1000 trial data sets and 
computed the fraction of simulated data sets in which the dif- 
ference of the SSFRs is larger than or equal to that based on 
the real data (fsim>obs)- The significance levels of the difference 
defined by 100(l-fsim>ob,s) (%)■ This test also confirms the in- 
creased SSFRs of interacting galaxies with late-type neighbors 
compared to non-interacting galaxies with a significance level 
of >99% for both strongly- and weakly-interacting galaxies in 
local and high-z universe. 

Similarly, the statistical significance of different SSFR dis- 
tributions of interacting galaxies with early- and late-type neigh- 
bors at R„ < 0.5rvir,nei IS examined with the K-S and MC tests. 
The difference is confirmed with a significance level of 98 and 
>99% from the K-S and MC tests, respectively, in the local uni- 



verse, and of 98% from both K-S and MC tests in high-z uni- 
verse. 

To study the evolution of SSFRs depending on the morphol- 
ogy of the nearest neighbor, we plot, in Fig.|5j the median SSFRs 
of strongly- and weakly-interacting galaxies as a function of red- 
shift in comparison with those for non-interacting galaxies. We 
plot the SSFRs of interacting galaxies normalized by those of 
non-interacting galaxies at each redshift bin in order to focus 
only on the relative ratio of SSFRs of interacting galaxies to 
those of non-interacting galaxies. 

For the early-type neighbor case (left panel in Fig. |5]l, 
the SSFRs of interacting galaxies are similar to or slightly 
smaller than those of non-interacting galaxies in all redshift bins. 
However, for the late-type neighbor case (right panel), it is seen 
that the increased SSFRs for weakly- and strongly-interacting 
galaxies are systematically larger than those for non-interacting 
galaxies in all redshift bins. The SSFRs of strongly- and weakly- 
interacting galaxies are, on average, larger than those of non- 
interacting galaxies by factors of about4.0±1.2 and 1.8+0.7, re- 
spectively. If these increased SSFRs are simply due to the blend- 
ing of two galaxies with similar masses in one FIR beam with- 
out any enhanced SFA, the SSFRs are expected to increase by 
a factor of only two. In fact, the neighbor-separation scale of 
Rn ~ 0.05rvii-,nei (~15 h^^ kpc for galaxies at Mr = -21) is 
important because the galaxies in pairs start to merge at this 
separation, so the SFA starts to change abruptly (Park & Choil 
2009). Thus the increased SFR for strongly-interacting galaxies 
(^n ^ 0.1rvii-,nei ~30 /i"' kpc for Mr = -21) by a factor of four, 
really reflects the enhancement of SFA due to a merger 

3.2. Change in dust temperature as a function of pair 
separation 

In Fig. |6] we show the dependence of dust temperature on the 
distance to and the morphology of the nearest neighbor galaxy. 
For local SDSS galaxies in (a), we plot the dust temperature 
derived from the flux density ratio (^eo/^ioo) using Eq. ([TJ to 
increase the statistics. For high-z GOODS galaxies in (b), we 
do not use the dust temperature derived from the FIR flux den- 
sity ratio because the large scatter in the correlation between 
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Herschel flux density ratios and dust temperature makes statis- 
tics worse. 

The dust temperature of both local and high-z late-type 
galaxies appears to increase as they approach late-type neigh- 
bors, which is clearly seen at /?„ S O.lrvir.nei- We also show 
the distribution of dust temperature in Fig.|4]for interacting and 
non-interacting galaxies in SDSS and GOODS. The MC test 
supports this different Tdust distribution between strongly- and 
non-interacting galaxies having late-type neighbors with a sig- 
nificance level of 96 and >99% for local and high-z galaxies, 
respectively. However, there are only two GOODS galaxies at 

< O.lrvir.nei, SO this trend needs to be checked with larger 
data set in future studies. On the other hand, the dust temper- 
ature of local late-type galaxies with early-type neighbors at 
Rn < O.lrvir.nei sccms to be lower than or marginally similar 
to that of non-interacting galaxies. Therefore, the difference in 
the dust temperature of strongly-interacting galaxies with late- 
and early-type neighbors seems to be significant, which is sup- 
ported by the MC test with a significance level of 98%. For high- 
z galaxies, there are few late-type galaxies having close early- 
type neighbors at R„ < O.lrvir.nei with Tdust measurements, so we 
can not study the trend. 



4. Discussion 

The increase of SSFR for galaxies in pairs found in this study is 
consistent with results in previous studies at low redshifts (e.g. . 



Kennicutt et a 



.1ll987[ iBarton et al. H2O00l . 
2004t IWoods & G elled 1200' 



Lambas et a\} '20031 



Nikolic et_alj l2004t |Woods & G elled 2007t iHuang & HwangI 
2011 ) and high redshifts (e.g.. lLin et al.l2007l:IWonget al.l20lir 
However, note that the amount of increased SSFR compared 
to isolated galaxies is different depending on the studies be- 
cause the definition of interacting galaxies and the observational 
selection effects are different. The increase of dust tempera- 
ture of local late-type galaxies strongly interacting with other 
late-type neighbors is also consistent with previous studies in 
the sense that the contribution of warm du st compared to cold 
dust increases with the merging sequence (iTelesco et al.lll988t 
[Xilouris et al. 2004). However, we found the hint of this trend for 
high-z galaxies at 0.4 ^ z ^ 1 .2 for the first time. It is also inter- 
esting to see that the dust temperature of local, late-type galaxies 
does not seem to increase when they are strongly interacting with 
early-type neighbors (see Fig. |6]l. Note also that the distance to 
the nearest neighbor might not be a direct measure of the merg- 
ing sequence because galaxies in pairs would merge after several 
encounters and its orbital geometry is complicated. 



We wish to emphasize that our analysis represents a major 
improvement compared to previous studies because we 1) have 
robust SFR measurements of galaxies with well-constrained FIR 
SEDs thanks to Herschel and AKARI, 2) have galaxy samples 
with high spectroscopic completeness (i.e. > 85% for SDSS 
galaxies at nir < 17.77 and 71 - 86% for GOODS galaxies at 
m- < 23), 3) study only the galaxies that are complete in terms 
of neighbor galaxies with stellar masses similar to target galax- 



ies (Mr.r 



A^r.tarset + 0.5), 4) distinguish the morphology of 



galaxies in pairs, which is an important indicator of the gas frac- 
tion in galaxies, and 5) use a distance to the nearest neighbor 
normalized by the virial radius of the neighbor instead of the 
physical distance in order to account for different masses and 
corresponding crossing time of galaxy samples. 

The dependence of SFR (or SSFR) on the morphology of the 
neighbor at R^ ~ O.Srvir.nei seen in Fig. |3]may imply that the hy- 
drodynamic interactions with the nearest neighbor play critical 
roles in triggering the SFA of galaxies in addition to the tidal in- 
tera ctions. This dependence is also seen for local normal galax- 
ies (iPark& Choi' '2009'; 'Xu et al.ll20To l) and local LlRGs and 
ULIRGs (Hwang et al. 2010^ and can be explained as follows. 
If a late-type galaxy approaches a late-type neighbor within half 
of the virial radius of the neighbor, the cold gas inflow into the 
central region of the target galaxy from the neighbor galaxy as 
well as from the disk of the target galaxy increases, which re- 
sults in the enhanced SFA of the target galaxy. Because of this 
starburst mode of compact star formation, the SSFR and the dust 
tempe rature ar e expected to increa se as we observed ( Elbaz et al] 
I2OIII; see also lChanial et al.ll2007h . Then when two galaxies fin- 
ish merging, the end product of the merger will be bright due 
to the very recent SFA, and the new nearest neighbor galaxy of 
the merger product will be far away. This may explain the exis- 
tence of some non-interacting galaxies at R^ > rvir.nei with large 
SSFRs and Tdust as seen in Figs.[3]and|6] 

On the other hand, if a galaxy approaches an early-type 
neighbor within the virial radius of the neighbor, the hot gas of 
the early-type neighbor prevents the galaxy from forming stars 
with cold gas, and/or there is no inflow of cold gas from its 
early-type neighbor, so the SFA of the galaxy is not boosted 
even if it has a close companion. The SF quenching mecha- 
nisms of hot gas in early-type neighbors could be similar to 
those of a hot intracluster medium of galaxy clusters acting on 
late-type galaxies in it, which are hydrodynamic processes such 
as thermal evaporati on, strangu lation, ram pressur e stripping, 
or viscous stripping (iPark & Cho i 2009; P ark & Hw ang 200^. 
Indeed, the X-ray observations of galaxies in pairs with mixed 
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morphology show evidence for extended X-ray halos of the early 
type that surround the late type, which supports our interpreta- 
tion (Grutzbauch et al. 2007). In addition, the fact that SSFRs in 
Fig. |5] are increased by similar factors in all redshift bins, sug- 
gests that similar physical mechanisms (such as hydrodynamic 
interactions plus tidal interactions described above) affect on the 
SPA over at least 8 billion years. 

The SFRs of galaxies are known to depend strongly on 
the local density, in the sense that spatially averaged SFRs of 
galaxies or star-forming galaxy fractions decrease as the back- 
ground densit^Jiicreas^^ the local u niverse (e.g., Lewis et al. 



200?, 'Gomez et all llool iParketalJ 120071: iLeeetalJ 120 lOt 
Hwang et al. 2010ah . In the high-z un iverse, some studies found 



similar results (e.g., iPatel et alj|2009l: iFeruglio et a"l]|2OI0h . but 
there are also hints for an opposite trend (i.e. increasing SFRs 
or star-forming galaxy fractions with increasing the background 
density) known as the reversal of SF R-density relation (e.g., 
lElbaz et al.l2007tlPopesso et alj20I lallbl) . When we consider the 
morphology-density relation ( Dresslei Il980l) . early(late)-type 
neighbors could be preferentially selected in high(low)-density 
regions, and the pair separation is correlated with the local den- 
sity. Therefore, the SFRs appear to depend on both large- and 
small-scale (attributed to the nearest neighbor) environments. 
One might expect that the difference in FIR properties depending 
on the morphology of and the distance to the nearest neighbor is 
simply due to the statistical correlation between the local den- 
sity and the properties of nearest neighbor. Thus, it is necessary 
to disentangle the effects of both environments. 

In previous studies, we have examined the effects of both 
large- and small-scale environments on the S FA and the mor- 
phologic al transfo rmation of SPSS ga laxies dPark e t al.*l2008t 
iPark & C hoi 2009: lHwang et al.ll20I0ah . and on the m orpholog- 
ical transformation of GOODS galaxies ("Hw ang & Parkll2009l) . 
We found that the SFA of local galaxies is still strongly affected 
by the nearest neighbors differently depending on the morphol- 
og y even at fixe d lar^e-scale backgr ound density (see F ig. 14 
in iHwang etal] |2010a and Fig. 7 in iPark & Choil |2009|) . The 
morphological transformation of high-z galaxies shows a sim- 
ilar trend (see Fig. 6 in Hwang & Park 2009). Therefore, it is ex- 
pected that the role of the large-scale environment in the SFA for 
high-z galaxies to be weak. We tried to investigate the change in 
SFRs of GOODS galaxies depending on both large- and small- 
scale environments, and found that galaxies with early (late)- 
type neighbors are not preferentially selected in high (low)- 
density regions. This implies that there is no bias in our results 
introduced by the difference in the large-scale environment (see 
also Fig. 6 in Hwang & Park 2009). However, we could not draw 
meaningful conclusions on the distinction of the effects of large- 
and small-scale environments on the SFA due to the small num- 
ber statistics, which needs to be investigated with larger data set 
of IR-detected galaxies in future studies. On the other hand, a 
detailed analysis focusing on the effects of large-scale environ- 
ment on the SFA of high-z galax ies can be found in other studies 
based on similar GOODS data (lElbaz et alJl2007t iPopesso et alJ 
l20TTillbl) . 



5. Conclusions 

Using the /fersc/ze/-selected galaxies in the GOODS fields and 
the IRAS plus AKARI-selected galaxies in the field of SDSS 
DR7, we studied the impact of galaxy-galaxy interactions on the 
FIR properties of galaxies and its evolution at < z < 1.2. Our 
main results follow. 



1. We find that the SFRs and SSFRs of galaxies, on average, 
depend on the morphology of and the distance to the nearest 
neighbor galaxy for all redshifts, within < z < 1.2. When 
a late-type galaxy has a close neighbor galaxy, the SFR and 
SSFR increase as it approaches a late-type neighbor, which 
is supported by K-S and MC tests with a significance level 
of >99%. However, the SFR and SSFR decrease or do not 
change much as it approaches an early-type neighbor The 
bifurcations of SFRs and SSFRs depending on the neigh- 
bor's morphology are seen at /?„ x O.Srvh^nei, which is also 
supported by K-S and MC tests with a significance level of 
>98%. 

2. For the redshift range < z < 1.2, the SSFRs of late-type 
galaxies having late-type neighbors are increased by factors 
of about 1.8+0.7 and 4.0±1.2, respectively, for the cases 
of weakly- (/?„ < O.Srvir.nei) and strongly-interacting (R„ < 
O.l^vir.nei) galaxics Compared to those of non-interacting 
galaxies. 

3. The dust temperature of both local and high-z late- type 
galaxies strongly interacting with late-type neighbors ap- 
pears to be higher than that of non-interacting galaxies with a 
significance leve of 96-99%. However, the dust temperature 
of local late-type galaxies strongly interacting with early- 
type neighbors seems to be lower than or similar to that of 
non-interacting galaxies. 

Our results suggest that galaxy-galaxy interactions and 
mergers have been strongly affecting the SFA and the dust prop- 
erties of star-forming galaxies over at least 8 billion years. 
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